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ABSTRACT 


The natural frequencies and mode shapes of a 1/15 scale 
space shuttle dynamics model are analytically determined. The 
model, a parallel beam type structure with delta wings, is 

— dynamically representative of -the— sti ffness -and -mass— proper±ies 

of an early space shuttle design. Important characteristics 

of the model are elastic interfaces with adjustable spring rates. 

Normal mode computations are made using the finite element 
modeling technique as implemented in the NASTRAN (NASA Structural 
Analysis) computer program. The feasibility of neglecting 
elastic deformations in the lower modes was investigated using 
a rigid body model. 

Using NASTRAN, natural frequencies and mode shapes were 
first calculated for the booster fuselage, orbiter fuselage, 
and both delta wings in a free-free condition. Next, the 
fuselages were connected for various spring rates. Then the 
wings were attached to each fuselage, and the booster and 
orbiter were analyzed as separate airplanes. Finally, the two 
airplanes were elastically joined, and the complete model 
was analyzed. 

In the complete model for an interface nominal spring 
constant of 10 3 lb/in there were forty modes with frequencies 
less than 200 HZ. The coupling of the booster and orbiter 
motions were found to be significantly affected by the spring 
interface stiffnesses. 
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NOMENCLATURE 


[a] 


height of orbiter attachment bracket, 
see Figure A-l. 

flexibility matrix. 


A 


cross-sectional area. 


A P 

A i 

C 

d 

E 

fi 

G 


cross-sectional area of one pitch leaf. 

coefficient in equation of constraint, 
see equation 4. 

width of booster attachment bracket, 
see Figure A-l. 

distance from top of booster attachment 
bracket to orbiter centerline. 

modulus of elasticity. 

forces applied at grid point i, see equation (1) 
modulus of rigidity. 


h 
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vertical distance between pitch leafs, 
see Figure A-l. 

area moments of inertia for CBAR elements. 
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area moments of inertia of booster attachment 
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area moments of inertia of pitch leafs. 

area moments of inertia of yaw leafs. 

mass moments of inertia of plane rigid body 
model, see equation (B-3) . 

mass moment of inertia of concentrated masses 
representing fuselage roll inertia. 

mass moments of inertia of rigid body model 
of orbiter. 

mass moments of inertia of rigid body model 
of booster. 
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polar moment of inertia, 
spring constant, 
stiffness matrix. 

length of pitch spring, see Figure A-l. 

total height of interface, see Figure A-4. 

length of yaw spring, see Figure A-l. 

effective lengths of booster attachment 
bracket arms for bending in XZ and YZ planes, 
respectively. 

concentrated mass. 

mass matrix. 

masses of plane rigid body model, see 
equation (B-3) * 

rigid body mass matrix of orbiter, see 
equation (B-8b) • 

rigid body mass matrix of booster, see 
equation (B-8c) . 

generalized displacement coordinates. 

rotational displacement components in mode 
shapes. 

rigid body transformation matrix, see 
equation (1). 

thickness of pitch leaf, see Figure A-l. 

thickness of yaw leaf, see Figure A-l. 

translational displacement components in 
mode shapes . 

displacement vector. 

strain energy of elastic interface. 

coordinate axes. 

\ 

coordinates of dependent grid point in 
equations of constraint, see equation (5) and 
Figure 7. 
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NASTRAN flexibility matrix. 

dimensions of wing attachment brackets, 
see Figure 4. 

transformation matrices, 
deflection components. 

rotation components about X-._o.r_ Y axes, 
mass density. 

rotation components about Z axis, 
circular frequency. 


NASTRAN notation: 


CBAR 

CELAS2 

C0NM2 


beam element connecting two grid points, 
scalar spring element connecting two grid points 
concentrated mass at a grid point. 


C0RD2R rectangular coordinate system defined by 

three points . 

GENEL general element connecting two grid points; 

defined in terms of flexibility matrix [Z] 
and transformation matrix [S] . 

MPC multipoint constraint; used to define rigid 

connections . 


0MIT1 

NSM 

PBAR 

SEQGP 


defines degrees of freedom to be deleted from 
problem through matrix partitioning; used to 
remove mass singularities. 

nonstructural mass. 

specifies properties of CBAR element. 

used to reidentify the sequence of grid points 
such as to optimize bandwidth of matrices. 
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INTRODUCTION 


The determination of the dynamic characteris tics of a 
spacecraft structure is an important step in the development of 
the overall spacecraft systdm. This is particularly true for 
the Space Shuttle since a new concept in space vehicle design is 
being used. The Space Shuttle will be comprised of a parallel 
arrangement of the launch vehicle and space vehicle during the 
launch phase. This arrangement raises the possibility of dynamic 
effects which have not previously been encountered. 

To initiate the dynamic studies of the Space Shuttle 
configuration, an unpublished analysis was made at LRC of a 
proposed design supplied by the NASA Manned Spacecraft Center. 

This analysis, for a straight-wing configuration of the Shuttle, 
represented the vehicles by an assembly of finite beam elements 
and lumped masses. Each of the two elastic connections between 
the vehicles was represented by a single translational spring. 

Pitch motions of the fuselages and flapping motions of the 
wings were permitted. Normal mode calculations were made for 
six different values of the connecting springs. The results 
of this ''exploratory study showed that, in general, the modes were 
significantly controlled by the spring rates of the vehicle 
connections. 

As the next phase in the study of the Space Shuttle dynamic 
behavior,. Langley Research Center designed and procured a 1/15 
scale model of proposed straight wing and delta wing configurations 
of the Shuttle. The vibration characteristics of this model are 
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being analytically and experimentally investigated. 

During an 11 week period at Langley Research Center in 
the summer of 1970, the principal investigator, under the NASA- 
ASEE Summer Faculty Fellowship program, was involved in 
computations to establish pre-test frequencies and mode shapes 
for the straight wing configuration of this model. These 
computations made use of the finite element modeling technique 
as implemented in the NASTRAN (NASA Structural Analysis) 
computer program. Some of the results of this work and some 
preliminary experimental results were presented at the NASA 
Space Shuttle Technology Conference held at LRC in March, 1971, 
(reference 1) . 

The present study is an extension of this work and is 
concerned with an analytical investigation of the delta wing 
configuration of this model. 

The objectives of this study were to: 

(a) Use the NASTRAN finite element computer program to compute 
the natural frequencies and mode shapes for the 1/15 scale 
dynamic model of the Space Shuttle with delta wings. 

(b) Study the effect of connection spring rates on frequencies, 
mode shapes, and the coupling between vehicle motion. 

(c) Formulate and analyze a rigid body representation of the 
1/15 scale model to determine the feasibility of neglecting 
elastic deformations in the calculation of the lower modes 
of the system. 
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DESCRIPTION OF MODEL 

In this section a brief description of the delta wing 
configuration of the 1/15 scale shuttle model will be given. 

The complete model will be- described first, and then the 
important characterist ics of each component will be reviewed. 

Complete Model 

The complete delta wing model is shown suspended by 
cables in Figure 1. This figure shows the booster component 
supporting the orbiter component through two elastic interfaces 
in a forward mounting position. It is also possible to support 
the orbiter through a rear mounting position which is not shown. 

A schematic with the orbiter mounted in the forward position 
is given in Figure 2. This figure also gives some of the 
general dimensions and defines the overall coordinate system 
used to describe the model. 

The stiffness and mass properties of the model are based 
upon the 12,500 lb. payload version of the MSC Space Shuttle 
Configuration. Each component was designed to closely approximate 
scaled stiffness and mass distribution curves of the MSC design. 

Fuselage Components 

The booster fuselage spar was fabricated from thin walled, 
tapered aluminum cylindrical sections reinforced by internal 
bulkheads. The orbiter fuselage spar is similar except that 
over a major portion of its length the diameter is constant. 

Each fuselage component has foam ballast weight holders used 
to support lead weights which simulate the propellant. The 
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orbiter fuselage includes a lead nose mass to simulate the 
payload. The locations of the weights are shown in Figure 2, 
and the values of the added weights are tabulated in Table I. 

Wing Components 

Each wing has three spars which were designed to give the 
proper bending stiffness distribution for the wing. The wings 
are clamped to the fuselages at three locations through 
mounting brackets welded to each spar. Over the length of 
each spar where the mounting bracket is attached the cross 
section of the spar is constant and has the form of an H. 

Slightly outboard of the mounting bracket the cross changes 
to a form obtained by welding two angle bars together. This 

cross section tapers linearly from the fuselage to the wing 
tip. The three spars are connected by forward-aft members of 

constant cross section. The forward-aft members are generally 
angle bars except at the wing tips where flat bars were used. 

Elastic Interface 

The details of one of the elastic interfaces are shown in 
a photograph in Figure 3. Each of the interfaces was designed 
to permit varying the nominal spring rates in the pitch and yaw 
directions. This was accomplished by incorporating removable 
pitch and yaw leaves. Four sets of spring leaves with nominal 
spring constants 10 2 , 10 3 , 10 4 , and 10 5 lb/in . were designed 
for each of the pitch and yaw directions. No scaling consider- 
ations were involved in selecting the nominal spring constant values. 
The values were selected from the results of preliminary calculations 
made at LRC. The dimensions of these leaves and other interface 
dimensions are given in Table A-I of Appendix A. 
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Fuselage-Wing Interface 

Some of the details of a typical booster and orbiter 
fuselage-wing interface are shown in Figure 4. The interface 
consists of two parts: (1) a top semi-circular bracket welded 
to the fuselage, and (2) a larger bottom bracket welded to 
the wing spar. The two brackets are joined by two bolts per 
side. The overall dimensions of the interfaces for the two 
fuselage are also shown in Figure 4. 

In this section a description of the model has been given 
In the next section the idealizations and properties of the 

mathematical model will be described. 
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ANALYSIS 


The analysis of the natural frequencies and mode shapes of 
the 1/15 scale model was made using a NASTRAN idealization and 
a rigid body idealization. 

In the NASTRAN analysis, each component was first analyzed 
in a free-free condition. As the second step, the fuselage 
components were connected by the elastic interfaces. As the 
third step, the wings were attached to each fuselage component, 
and the booster and orbiter were analyzed as separate airplanes. 
In the final step, the two airplanes were connected through the 
elastic interface to represent the complete model. 

In the rigid body idealization the plane motion of the two 
elastically connected fuselages was analyzed first. Then a 
three-dimensional analysis of the complete model was made. 

NASTRAN Analysis 

In this section the NASTRAN idealization of the model 
will be described. The idealization of each component will 
be described first. Then two methods of representing the 
elastic interfaces will be described. Finally, a description 
of the method of representing the fuselage-wing interfaces 
will be discussed. 

Fuselage Components .- The booster fuselage was idealized by 
an assembly of ten concentrated masses and nine beam elements. 
The concentrated masses were represented in NASTRAN using the 
bulk data card C0NM2 , and the beam elements were represented 
by CBAR elements. The idealization and the associated mass and 
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element properties of the booster fuselage are given in Tables 
II and III. 

In NASTRAN, the beam element CBAR is assumed to have uniform 
properties. Because of this the cross-sectional properties 
tabulated in Table II are average values based on the diameter 
at the mid-point of each element. 

Added weights 2 and 4 (See Figure 2) being distributed 
over relatively short lengths, were represented as concentrated 
masses at grid points 5 and 8, respectively. Added weights 1 
and 3 were distributed uniformly over elements 2 and 3 and 
are shown as non-structural mass (NSM) in Table II. (Non- 
structural mass is a quantity of material which adds mass to a 
structural element but does not contribute to the stiffness.) 

The actual mass of each beam element is computed internally in 
the NASTRAN program from the input data shown in Table II. 

The NSM property is input on the PBAR beam property card. 

Concentrated masses were also used to model local 
distributions of mass, such as the reinforcing bulkheads. The 
values of mass and mass moments of inertia shown in Table III 
for the Liftoff condition are listed as nominal values and 
test values. The nominal values were computed from the 
drawings of the model. The test values were obtained by 
adding laboratory weighed values of the elastic interfaces 
and rigging brackets to the nominal values. 

The mass moment of inertia I xx about the roll axis of the 
fuselage was represented by concentrations at each grid point. 
Included in each value of I xx are the roll inertias of the 
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bulkheads, the added weights, and each element of the fuselage. 
The inertia of each element was divided equally between adjacent 
grid points. In Table III for Liftoff, two values of the roll 
moments of inertia 1^ are shown. The nominal values are 
based upon the design dimensions of the model. The test values 
were calculated for slightly different radial locations of 
the weights used in the experimental vibration tests. The 
mass moments of inertia of the add-on weights and the fuselage 
itself about the other two axes, i.e., the "rotatory inertia", 
was neglected. 

The orbiter fuselage was idealized in a manner similar to 
the booster fuselage. Seven beam elements and eight 
concentrated masses were used. The idealization and properties 
of the orbiter fuselage are shown in Tables IV and V. 

The two added weights of this component were distributed 
over their supporting beam elements and modeled as NSM. The 
payload was modeled as a concentrated mass. 

For all elements a modulus of elasticity E =10.1 x 10 6 
lb/in 2 and a modulus of rigidity G = 3.8 x 10 6 lb/in 2 were 
used. The mass density was taken as p = 2.539 x 10 -4 lb-sec 2 /in 4 . 
These values were also used for the wing components. 

Wing Components .- The booster delta wing is shown 
schematically in Figure 5. Three coordinate systems were used 

to locate the grid points and are shown in this figure. The 

coordinates of the grid points in these coordinate systems 

t ' 

are tabulated in Table VI. In NASTRAN, coordinate system 1 was 
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used as the basic coordinate system, and coordinate systems 
2 and 3 were referenced to the basic coordinate system by 
using the bulk data cards C0RD2R. 

The wing was idealized using 55 CBAR elements. For the 
spars which are tapered the average properties of the elements 
were used. The properties for all elements are tabulated in 
Table VII. All properties listed in this table are referred 
to the local beam element coordinate system. II represents 
the moment of inertia of the cross section associated with 
bending in a plane perpendicular to the plane of the wing; 

12 represents the moment of inertia associated with bending 
in the plane of the wing. 112 represents the product of 
inertia in the local NASTRAN coordinate system. The sign was 
established according to the NASTRAN sign conventions given 
in paragraph 1.3.2 of Reference 2. 

The lower half of the fuselage mounting brackets which 
are welded to the wing spars were represented as concentrated 
masses. Equal masses of 10.34-4 lb-sec 2 /in were located at 
grid points 1-9. The bending stiffnesses of these brackets 
were neglected. 

The orbiter delta wing is shown schematically in Figure 6. 
The grid points for this wing were located in one basic 
coordinate system, and these coordinates are tabulated in 
Table VIII. 

The wing was idealized using 29 CBAR elements. As for the 
booster wing, average values were used for the tapered portions 
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The element properties are tabulated in Table IX. The definitions 
for these quantities are similar to those for the booster wing. 

Equal concentrated masses of 4.66-4 lb-sec 2 /in were located 

at grid points 1-9 to represent the fuselage mounting brackets. 

Elastic Interface .- The' elastic interfaces were modeled in 

two ways. In the first model the interface was represented using 

simple linear springs with the nominal spring constants in the 

pitch and yaw planes. In NASTRAN these are called scalar elements 

and were input as bulk data through the use of the CELAS2 card. 

In the second representation of the elastic interfaces the 

general structural element GENEL was used. The element describes 

the elastic interface through a six by six flexibility matrix [Z] 

and a six by six rigid body transformation matrix [S] . These 

matrices are defined by the equation, (See Reference 3, Section 5.7), 

{ uj_ } = [Z] {fi> + [S]{u d > (1) 

where {u^} denotes the set of displacement components 

at grid point i; 

r Z ] denotes the flexibility matrix; 

{f^} denotes the forces at grid point i corresponding 
to the u^ displacements, 

[S] a rigid body matrix whose terms are the displacements 
(ui) due to unit motions of the coordinates (u d ), 
when all f ^ = 0 . 

{u^} denotes the set of displacement components at 
grid point j. 

For example, the forward elastic interface connects grid 
point 3 on the booster fuselage and grid point 22 on the orbiter 
fuselage. In this case 

(ui) = [22- 1, 22-2, 22 - 3 , 22 - 4 , 22-5, 22-6 J , T 

and i T ( 2 ) 

{u d } = [_3-l, 3-2, 3-3, 3-4, 3-5, 3-6J, 
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where 22-1 denotes the displacement component 1 at grid point 
22, 22-2 denotes the displacement component 2 at grid point 22, 
etc. ' J 

With no applied forces {f-p} the displacements (u^) and {ud) 
are related by the transformation matrix [S] . For the geometry 
show n, in F i g ure 2 , 


[SJ = 


1 

0 

0 

0 

0 

0 


0 

1 

0 

0 

0 

0 


0 

0 

1 

0 

0 

0 


0 

■15.6 

0 

1 

0 

0 


15. 6 
0 
0 
0 
1 
0 


0 

0 

0 

0 

0 

1 


(3) 


The flexibility matrix [ Z ] was determined for each of the 
nominal spring constants using the theory developed in Appendix A . 

The results of these calculations are tabulated there in Table A-II. 

In addition, for the nominal spring case of 10 3 lb/in, the 
flexibility coefficients were measured experimentally at Langley 
Research Center. The non-zero elements of the theoretical 
flexibility matrix for the nominal 10 3 lb/in case and the 
experimental values for this case are tabulated in Table X. 

Fuselage-Wing Interfaces .- Each fuselage-wing interface was 
assumed to be rigid. In NASTRAN rigid connections are modeled 
by writing equations of constraint between the displacement components 
at the grid points of interest. These equations take the form 


l A j u j= 0 


(4) 


where Aj are coefficients determined from geometry and Uj are 

\ 
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the displacement components. The first coordinate in the sequence 
is assumed to be the dependent coordinate. Equations of the form 
of (4) are called multipoint constraints and are input as bulk 
data using MPC cards. 

For the fuselage-wing interfaces a general set of multi- 
point constraint equations were derived as follows. Consider 
grid points A and B shown in Figure 7. For convenience, the 
origin of the coordinate system is located at grid point A. 

Grid point A and B are connected by a fictitious rigid link. The dis- 
placements at point B are then related to the displacements at 
A by the equation 



Equation (5) represents 6 equations of constraint each of the 
form of equation (4). The multipoint equations of constraint 
for NASTRAN were written by substituting appropriate values of 
x 0 , y 0 , z 0 in these equations. 

The wings were attached to the fuselage by assuming a rigid 
connection between six grid points at the outboard edge of the 
mounting brackets on the wing spars and the nearest grid points 
on the fuselage. Thus the wings were "cantilevered" from the 
fuselage with the cantilever points below and outboard of the 


fuselage centerline. 






For the booster fuselage-wing interface grid points 6, 8, 

9 on the booster fuselage were connected to grid points 34-39 
on the booster wing. These connections are shown schematically 
in Figure 8. The values of x 0 , y 0 , z 0 used in equations (5) are 
tabulated in this figure. A total of 36 equations of constraint 
were written for the Booster fuselage-wing interface. 

The orbiter fuselage-wing interface was modeled in a 
similar way and is shown schematically in Figure 9. For this 
interface, 36 equations of constraint were also written. 

Booster and Orbiter Airplanes .- The booster fuselage and 
delta wing components were connected using the multipoint 
equations of constraint to give a complete airplane. The same 
grid point numbers were used for the fuselage; new grid points 
for the wings were obtained by adding a double zero suffix. A 
plan view of the booster airplane idealization is shown 
schematically in Figure 10. The interface masses were equally 
distributed at grid points 34-39 of this idealization. 

The orbiter fuselage and delta wing components were 
connected in a similar way to yield a complete airplane as shown 
in Figure 11. New wing grid point numbers were obtained by adding 
a zero one suffix to the original numbers. The interface masses 
were equally distributed at grid points 44-49 of this model. 

Complete Model .- Using two GENEL elements, the booster and 
orbiter airplanes were elastically connected ito represent the 
complete model. The forward GENEL element connected grid point 
3 on the booster fuselage to grid point 22 on the orbiter fuselage 
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The aft GENEL element connected grid point 4 on the booster 
fuselage to grid point 26 on the orbiter fuselage. 

To reduce the band width of the stiffness matrix for the 
complete model the grid points were renumbered. This grid point 
resequencing was input to NASTRAN through the SEQGP bulk data 
cards. The fesequencing of grid points is tabulated in Table XI. 

Rigid Body Analysis. 

Two mathematical idealizations were used for the rigid 
body vibration analysis. In the first idealization the motion 
was assumed to be restrained to the pitch plane with longitudinal 
displacements prohibited. This gave a simple four degree of 
freedom model. The elastic interface for this case was 
represented as a scalar spring. In the second idealization, 
general three dimensional motion was permitted. The booster and 
orbiter were each allowed to have three displacements and three 
rotations so that the system possessed twelve degrees of freedom. 
The elastic interfaces in this case were represented by twelve 
by twelve stiffness matrices. 

The details of these analyses are presented in Appendix B. 
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RESULTS AND DISCUSSION 


In this section the results of the analysis will be 
presented and discussed beginning first with frequencies and 
mode shapes of each component. Then the results for the 
elastically connected fuselages will be described. Next the 
modes of the booster and orbiter airplanes will be discussed, 
and finally, the frequencies and mode shapes for the complete 
model will be given. 

Component Modes 

Fuselages .- The natural frequencies of the booster and 
orbiter fuselages were calculated using the properties given 
in Tables II-V for the nominal masses and moments of inertia. 

The booster fuselage has 60 degrees of freedom (DOF) with all 
displacements permitted; the orbiter has 54 DOF permitting all 
displacements. The consistent mass matrix option was specified, 
and the eigenvalue problems were solved by the inverse power 
method. The natural frequencies of the booster and orbiter 
fuselages are tabulated in Table XII. For comparison, some pre- 
liminary experimental frequencies obtained at LRC are also listed 

For the nominal weights these results show for the booster 

there are one to three bending modes in the frequency range 

below 200 HZ depending on the weight condition. For the cases 

where the weights were added, a torsional mode was also 
predicted at about 114 HZ. The axial modes were found to be 

outside of the frequency range of interest. For the orbiter 

fuselage, the frequencies tended to be substantially higher 

than the booster fuselage with only the fundamental bending 


mode below 200 Hz. 
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The experimental frequencies tend to generally confirm 
the predicted values particularly in the lower modes. The 
agreement becomes poorer with an increase in modal number. 

Better agreement with the experimental values for the higher 
modes could possibly be found by including the rotatory inertia 
of the added weights and by further subdivision of the fuselages 
using more beam elements. Because of the reasonable agreement 
of the lower modes these steps were not pursued. 

The mode shapes of the bending modes of each fuselage are 
shown in Figures 12 and 13 for the liftoff and burnout weight 
distributions. These modes correspond to either pitch or yaw 
bending due to symmetry. 

Wings . - In the "free-free" condition with all grid points 
having 6 displacement components the booster delta wing has 198 
DOF, and the orbiter delta has 162 DOF. To reduce the size 
of the eigenvalue problems the lumped mass option was selected 

•<r 

for each wing, and the rotational displacement components were 
deleted from the eigenvalue problem through the use of the bulk 
data card 0MIT1. This reduced the size of the eigenvalue problems 
for the booster and orbiter wings to 99 and 81 DOF, respectively. 
Both of these eigenvalue problems were solved using the Givens 
method. The natural frequencies for the two wing components 
are listed in Table XIII. Some preliminary experimental 
frequencies are also listed. 

The frequency tabulations show that the wings have relatively 
dense frequency spectrums compared to the fuselage components. 

The booster wing has 20 frequencies below 200 HZ, and the orbiter 

wing has 16 frequencies in this range. 
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The experimental frequencies show reasonable agreement with 
the predicted values; however, a large number of the theoretical 
modes were not detected in the vibration tests. This may be 
due to a number of reasons including the fact that a number of 
the frequencies tend to be very closely spaced. Also, it is 

likely that the higher modes of the analytical^model are„ some- 

what in error because of the crudeness of the finite element 
model. 

There are at least two ways in which the mathematical model 
could have been improved. In the central region of eadh wing 
the mounting brackets are welded to the spars. Although the 
mass of these brackets was included, the stiffening effect of 
these brackets was neglected. No further attention was given 
to improving the finite element model of the wing in this region, 
however, because for the assembled components the fuselage-wing 
interface was assumed to be rigid. The second way of improving 
the finite element model concerns the tapered spars. Outboard 
of the fuselage brackets the wing spars taper to give a 
relatively rapid variation of El. For simplicity, one finite 
element of constant cross section was used to represent each 
spar between junctions with the fore-aft wing members. For 
better agreement in the higher modes these sections of the 
wing should be represented with more elements. Since this 
would have roughly doubled the DOF for each wing further sub- 
division of the elements was not attempted. 
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. The NASTRAN structural plotting capability was used to 
generate plots for the wing components. Figure 14 shows the 
plotter view of a typical delta wing. In the plots, each wing is 
viewed from an aft position at an angle of about 23° above the 
X axis. The first nine modes of the booster wing viewed in 
this manner are shown in Figure 15; the first nine modes of 
the orbiter wing are shown in Figure 16. 

Each wing is symmetrical about the centerline of the 
fuselage, and the modes are therefore either symmetrical or 
antisymmetrical about this axis. Most of the modes are 
characterized by flapping motions, i.e. displacements which 
are primarily perpendicular to the plane of the wing. For 
the booster wing, of the mode shapes shown, all modes except 

3 and 5 illustrate this type of deformation. Modes 3 and 5 
consist of fore and aft deformations in the plane of the 
wing. For the orbiter wing, all modes shown, except for modes 

4 and 7, have primarily out-of -plane displacements. Mode 4 
has mainly in-plane displacements, and mode 7 has a large in- 
plane displacement at the centerline but the wing tips move 
perpendicular to the plane. 

Although most wing modes are characterized by a large 
displacement component in one direction, the modes are generally 
three dimensional. For this reason, it should be observed that 
plots of the mode shapes from a single view such as used in 
Figures 15 and 16 do not show the general complexity of the 
motion. 
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Elastically Connected Fuselages 
Nominal springs .- The frequencies of the elastically 
connected fuselages were first calculated with NASTRAN assuming 
plane motion with axial displacements prohibited. A consistent 
mass matrix was specified, and the inverse power method was 

used to solve the- eigenvalue problem,. The frequencies of the 

first six modes for the nominal weight condition and the four 
nominal spring cases are tabulated in Table XIV. 

Table XIV Frequencies of elastically connected fuselages. 

Nominal springs. Liftoff weight distribution. 


Frequencies, HZ 



K 

lb/ in 



Mode 

10 2 

10 3 . 

. 

s 

10 5 


3.53 

10.96 

27.50 

32.21 


10.16 

30.15 

41.24 

57.60 


36.90 

38.61 

77.89 

117.1 

■. ' ■ 

103.9 

107.3 

117.1 

156.0 

■ 

117.3 

117.4 

147.3 

196.9 

1 

198.8 

199.0 

201.5 

223.5 


The mode shapes for the nominal 10 2 lb/in spring are shown 
in Figure 17, and the mode shapes for the 10 s lb/in in spring are 
shown in Figure 18. 

For the 10 2 lb/in spring. Figure 17 shows the fuselages 
move as rigid bodies in modes 1 and 2. Moreover , modes 3,4, 
and 5 are essentially uncoupled fuselage bending modes for this 
spring constant. However, as the nominal spring constant is 
increased the fuselages have coupled deformations in all modes. 
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, This is clearly shown in Figure 18. 

The plane motion rigid body analysis was used to predict 
the frequencies for these spring and weight conditions, and the 
results are compared with the NASTRAN frequencies in Figure 19. 
This figure shows the rigid body analysis gives good agreement 
with the NASTRAN results for the first two modes for spring 
constants less than 10 3 lb/in. For spring constants greater 
than 10 3 lb/in., the rigid body analysis diverges rapidly 
from the NASTRAN results. Thus for K > 10 3 lb/in., the 
flexibility of the fuselage components play an important role 
and cannot be neglected. The mode shapes for the 10 2 lb/in 
case are compared in Figure 20. 

Interface matrix representation .- The frequencies for the 
elastically connected fuselages were also calculated with 
NASTRAN by representing the elastic interface by its theoretical 
flexibility matrix. The lumped mass approach was taken and 
rotational components of displacement in the XZ and YZ 
planes were omitted. Each grid point was thus permitted to 
have all three translational components of displacement and a 
rotational component of displacement about the X axis. The 
resulting eigenvalue problem involved 72 DOF and was solved 
by the Givens method. The test weights were used for these 
calculations. The frequencies for the first fourteen modes for 
the four nominal spring cases are tabulated in Table XV. 

The predominate motion associated with each frequency is 
also specified in this table. 
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With three dimensional motions permitted, the modes become 
more complex. Generally, the inodes may be identified with 
motion in either the XZ (pitch) plane or YZ (yaw) plane. However, 
as Table XV shows, combined pitch-axial modes and roll modes 
also exist. 

Computer plots of these modes were made for the 10 3 
lb/in case. The plotter view of the elastically connected 
fuselages is shown in Figure 21a. Significant displacement 
components of each end of the fuselages are given in each plot. 

The positive conventions of the displacement components are 
shown in Figure 21b. The first nine modes for the 10 3 lb/in 
case are plotted in Figure 22. 

Modes 1-6, 8, 9, and 11 of Figure 22 have essentially 
uncoupled motion in either the yaw or pitch planes. Mode 7, 
however, illustrates a strong coupling between axial and 
pitch plane displacements while mode 10 is characterized by 
coupling between yaw displacements and roll of the orbiter 
fuselage. 

For the nominal 10 3 lb/in case, the frequencies were also 
calculated using the experimentally determined flexibility 
matrix. These results are compared with the frequencies for 
the nominal spring constants and the theoretical flexibility 
matrix in Table XVI. The agreement between the frequencies 
predicted by the two different flexibility matrices was reasonably 
close. The only significant difference occurs for the frequency 
of the combined pitch-axial mode which the theoretical interface 
matrix predicts to be 75.2 HZ, and the experimental interface 
matrix predicts as 64.3 Hz. 
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As Table XVI indicates , the nominal interface predicts 
the frequencies of the pitch and yaw modes to be identical. 

The pitch frequencies for the first four modes were predicted 
fairly well, but the yaw modes and the combined pitch-axial 
modes were poorly predicted'. 

Analysis of the elastically connected fuselages shows 
the connection of the two fuselages produces significantly 
more frequencies with more complex mode shapes than the 
individual components possess. This is illustrated in a 
frequency spectrum plot in Figure 23a where the number of 
frequencies for the connected fuselages and the fuselage 
components are presented. The connected fuselages have 
thirteen frequencies less than 150 HZ while the individual 
components possess only four frequencies in this range. 

Some preliminary experimental data obtained at LRC are 
plotted in Figure 23b for comparison with the predicted 
values. The experimental data also shows rather dense frequency 
spectra at two widely separated frequency bands. Although some 
experimental and theoretical frequencies appear to be in agreement, a 
comparison of mode shapes will be necessary before further conclusions 
can be reached concerning the theoretical-experimental correlation. 

The first six modes of the elastically connected fuselages 
were also calculated using the rigid body analysis and a 
stiffness matrix representation of the interface. Some of 
these results are compared with NASTRAN calculated frequencies 
in Figure 24. Up to K = 10 3 lb/in, the rigid body analysis 
predicts the frequencies of the first four modes reasonably well. 
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Frequencies of the fifth and sixth modes, however, are 
predicted with considerable error and are not shown. Above 
10 3 lb/in the rigid body analysis results in considerable 
error even for the first or second modes. These observations 
are in agreement with the mode shapes shown in Figure 22 where 
it may be seen^that the fuselages undergo significant deformation 
after the fourth mode. 

Booster and Orbiter Airplanes 

The booster and orbiter airplanes were both analyzed using 
the lumped mass approach. To avoid a singular mass matrix all 
rotational coordinates on the wings were omitted on both 
airplanes. On the fuselages, the rotations in the XZ and the 
YZ planes were omitted at all grid points except the wing 
attachment points. Initially it appeared necessary that 
rotational components in the XZ and YZ planes at the wing attach- 
ment points should also be omitted. This was attempted and 
resulted in an ill-conditioned mass matrix. Subsequent 
investigation showed that the multipoint constraint equations 
removed the apparent mass singularity, and the rotational 
components at these points should not be omitted. Roll rotations 
were permitted on both fuselages. The booster airplane was 
described by 136 DOF, and the orbiter airplane had 110 DOF. 

Both eigenvalue problems were solved by the Givens method. 

The first thirteen frequencies of the booster airplane and 
orbiter airplane are tabulated in Table XVII. Computer plots 
of the mode shapes were used to classify the modes according 
to two types of predominant motion. The plotter orientation 
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for the airplane components and complete model is shown in 
Figure 25. Plots of six booster airplane modes are shown 
in Figure 26, and similar plots for the orbiter are shown 
in Figure 27. 

Both the booster and orbiter were characterized by the fact 
that only a few modes represented coupling between fuselage 
and wing displacements. The majority of the modes consisted 
only of wing motions with little or no fuselage deformation. 

The wing-type modes were either symmetric or anti-symmetric 
about the XZ plane. 

For the booster, modes 1 and 2 in Figure 26 illustrate 
coupled motions between the fuselage and wings. The next four 
modes shown consist solely of either symmetric or anti-symmetric 
wing deformations. Similar behavior for the orbiter modes is 
shown in Figure 27. In this case, however, the first four 
modes consist of only wing deformations and the fifth and 
sixth modes show fuselage-wing coupling. 

Complete Model 

The final representation of the complete model was obtained 
by joining the two airplane models with the matrix representation 
for the elastic interfaces. The lumped mass approach was taken 
with the same rotations deleted as in the airplane analyses. 

Using the Givens method a matrix eigenvalue problem with 246 DOF 
was solved. The total Computational time (CPU) on a LRC CDC-66Q0 
computer was 2900 seconds. 

The frequencies for the nominal 10 3 Ib/in interface are 
tabulated in Table XVIII. There were forty modes less than 200 HZ 
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and the frequencies of thirty-two of these are tabulated. 
Computer plots of the mode shapes were used to classify the 
modes according to their predominant motion. The plotter view 
of the undeformed complete model is shown in Figure 25b, and 
plots of ten selected mode shapes are shown in Figure 28. 

Of the thirty two modes tabulated, bight modes involve- 

coupled motions between the fuselages and twenty-four modes 
involve uncoupled wing motions. Mode 2 (Figure 28b) is a good 
example of a coupled fuselage-wing mode. The uncoupled wing 
modes occur with either the booster wing vibrating, the 
orbiter wing vibrating, or both wings vibrating simultaneously. 
Mode 13 (Figure 28h) shows the booster wing vibrating 
separately, mode 5 (Figure 28e) shows both orbiter wing 
vibrating separately, and mode 6 (Figure 28f) shows both wings 
vibrating in phase in a symmetric mode. Wing modes also occur 
in which one wing is vibrating symmetrically, and the other is 
vibrating anti-symmetrically . These are designated as combined 
wing modes in Table XVIII. Mode 14 (Figure 28i) is an example 
of this type mode. 

The complete model was also analyzed for the nominal 10 3 
lb/in case using the experimental flexibility matrix to 
represent the interface. The agreement between the results was 
good with all frequencies and mode shapes matching very closely 
except for one mode. The experimental flexibility matrix 
introduced a mode at 59.7 HZ not predicted by the theoretical 
interface. The mode consisted of coupled fuselage-wing 
deformations. 


25 



The connection of the booster and orbiter airplanes gives 
a relatively dense frequency spectrum. The theoretical 
frequency spectrums for the two airplanes and the complete model 
are compared in Figure 29a-29 c. 

Connecting the two airplanes introduces a number of new 
modes into the frequency spectrum. Examples are the coupled 
fuselage modes which occur at about 11, 14, 21 and 24 HZ. 
However, a "number of modes in the complete model may be 
identified with matching modes in the airplane components. For 
example, the booster wing mode at 65.6 HZ (Figure 28h) for the 
complete model can be identified as the booster airplane wing 
mode at 65.2 HZ shown in Figure 26f. A number of other modes 
also can be matched. However, it should be observed that 
earlier calculations for the fuselage components have shown 
that the spring interface stiffness plays an important role in 
the coupling of motions between individual components. Thus 
it can be reasonably expected that higher spring rates would 
give a considerably larger proportion of modes with coupled 
motions between the two components . 

Some preliminary vibration test data for the complete model 
are shown in Figure 29d. As with the fuselage components a 
number of frequencies tend to be in agreement, but a comparison 
of mode shapes will be necessary to establish definite 
correlations. 

The rigid body analysis was also performed for the complete 
model. For the 10 3 lb/in 2 interface the rigid body analysis 
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predicted the first mode frequency to be 12.1 HZ which is 11% 
higher than the NASTRAN value. The frequencies of the higher 
modes were in considerable error. This result was to be 
expected since Figure 28 from the NASTRAN analysis showed 
significant deformations in all modes. 



CONCLUSIONS 


On the basis of the vibration analysis of the 1/15 scale 
dynamic model of the Space Shuttle with delta wings the 
following conclusions were reached: 

1. The model has a relatively dense frequency spectrum 
with about forty modes less than 200 HZ. The mode shapes are 
generally three-dimensional with coupled fuselage-wing 
deformations. However, numerous uncoupled wing modes were 
present. 

2. The coupling of the booster and orbiter motions is 
affected very significantly by the spring interface stiffnesses. 
A good analytical representation of the interface is necessary. 
The representation of the interface by its flexibility matrix 
through the NASTRAN general element appears to be a satisfactory 
approach. 

3. The correlation of the predicted frequencies and 
preliminary LRC experimental frequency data was generally fair. 

A final conclusion on this correlation can be reached only 
after detailed comparison of the mode shapes. Some regions of 
the analytical model may require further refinements in the 
finite element representation. 

4. It is feasible to use a rigid body analysis to predict 
the lower modes only for very weak spring interfaces. Generally, 
for a nominal spring of less than 10 3 lb/in., a rigid body 
analysis is satisfactory for predicting only the first mode. 
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APPENDIX A 


Calculation of Elastic Interface 
Flexibility and Stiffness Matrices 

In this appendix the analytical determination of the 
flexibility and stiffness matrices of the elastic interface 
will be described. The elements of the flexibility matrix were 
determined for use in the general structural element GENEL in 
the NASTRAN program. The stiffness matrices were developed for 
use in the FORTRAN program for the rigid body vibration analysis. 
In the first part of this appendix the determination of the 
flexibility matrices will be described. Next the form of the 
stiffness matrix will be described in terms of its submatrices. 
Finally, relations between the stiffness and flexibility 
coefficients will be derived. 

Flexibility Matrix 

A schematic of the elastic interface is shown in Figure A-l, 
and the positive conventions for the displacements and rotations 
are shown in Figure A-2 . 

The flexibility matrix will be defined in terms of the 
displacements and rotations at the orbiter centerline. For 
calculating the flexibility matrix, the interface will be assumed 
rigidly clamped at the booster fuselage. The flexibility matrix 
then will take the form 
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u 5 u 3 u 9 u 7 U 1I 





U 1 

| a 55 a 53 | 

1 


^5 

j a 3 5 a 3 3 

1 


u 3 


r 

a 99 

n 

a 97 , 

u 9 

1 

I 

a 79 

a 77 |^ j 

u 7 



1 a T 

1 J 

U 11 


(AI-1) 


where all terms not shown are zero. The flexibility coefficients 
may be identified as follows: 

The coefficient a z may be identified with the 
axial displacment of the interface, i.e., the 
displacement in the Z direction; 

The two by two matrix with subscripts 5 and 3 is 
associated with displacements in the XZ plane; 

The two by two matrix with subscripts 9 and 7 is , 
associated with displacements in the YZ plane, and 
The coefficient a T is associated with torsional 
deformation of the interface, i.e., with the 
rotation about the Z axis. 

Equations for each of these coefficients will now be derived. 


Axial Displacement: 

The coefficient a* may calculated by applying a unit 

u 

force to the interface at the orbiter attachment point and 
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computing the displacement U^ (See Figures A 1 and A 2 ) . For this 
calculation it is reasonable to assume the axial displacement 
of the yaw spring is negligible as well as the axial displace- 
ments of the booster and orbiter brackets. Then the only 
displacement contributing to a z is due to the bending of the four 
pitch leaves in the YZ plane. Thus 


& 3 P 

4 8EIpx 


(A- 2 ) 


where I px represents the moment of inertia of one pitch leaf 
about an X axis through its centroid. 

Displacements in the XZ Plane: 

The flexibility coefficients 355 and 835 are found by 
applying a unit force at the orbiter centerline in the X direction 
and calculating the U5 displacement and the U3 rotation. 

In calculating the 115 displacement the following contributions 
were considered: 


( 1 ) Bending deflection of the yaw leaves as cantilever beams, 

(a+h)£y 


6 * " ZeI 


0 3 
x 


+ sm 


(A- 3 a) 


L YWY ‘* C '- L YWY 

where Iy^ represents the moment of inertia of one yaw leaf 
about a Y axis through its centroid. 

( 2 ) Deflection at orbiter cdnterline due to end rotation of 
the yaw leaves as canti levers, 

7 1 ~ 1 u 1 2 t 

t-Y 


(a+h)l z Y (a+h) 2 ^ 

- a ®. - - — : — - ■* 


(A- 3 b) 


4EIywy 2EIywY 

where 0 a is the end rotation at the yaw leaves as cantilevers. 
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(3) Deflection at orbiter centerline due to rotation of 
pitch leaves about a local Y axis, midway between the upper and 
lower pitch leaves. 


h 2 

. ■, (a+ j ) l 3 p 

<S 3 = aQ, = -A 2 p 


12 


h 2 


El 


(A-3c) 


PZ 


where Ip Z represents the moment of inertia of one pitch leaf 
about a Z axis through its centroid. — - 

(4) Bending deflection of booster bracket. For this 
calculation ’ the bracket was assumed to act as a cantilever 
beam with an effective length for bending in XZ plane. 


«. * 


_V + Jhl 

6 EI by 4ei BY 


(A-3d) 


where Igy denotes the moment of inertia of one arm of the 
bracket about a Y axis through its centroid. The quantity d 
is the distance from the top of the bracket to the orbiter 
centerline. 

(5) Deflection at orbiter centerline due to end slope of 
booster bracket. 




d 0 3 


+ 

4EI B y 


d * L i 

2EI B y 


(A-3e) 


The total deflection due to a unit force is 

a 55 = ^ i + ^ 2 + *5 3 + <5 4 + 6 5 . (A-3f) 

The rotation U 3 due to this unit load is assumed to be 
made up of the following components: 

(1) The end rotation of the yaw springs as cantilever beams, 

(A-4a) 


G i " 4EI 


Y ( a+h) Zy 

+ 2S1ywy 


YWY 
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(2) The rotation of the pitch springs '.about the local Y axis, 

{A- 4b) 


0 . 


/ _ . h. £ 3 
= ( a+y) p 


12h 2 EI pz 


(3) The end rotation of the bracket arms as cantilever beams, 


4EI 


Lj 2 dL x 


BY 


2EI 


(A-4c) 


BY 


The total rotation due to a unit force is 

‘ a 35 = 0 X + 0 2 + 0 3 • (A-4d) 

The flexibility coefficients a^ and a^ are found by 
applying a unit moment about a Y axis through the orbiter 
attachment point. By symmetry, 

a 53 = a 35 * (A-5) 


The rotation u 3 due to a unit moment has three components 

similar to the rotation components considered in a 35 . Thus 

Lx 


a 33 


2EI 


Ax. 

YWY 




12h z EI pz 2 E! by 


(A-6) 


Displacements in the YZ plane: 

The flexibility coefficients agg and a-/$ are found by 
applying a unit force in the Y direction and calculating the ug 
displacement and the u^ rotation. 

In calculating the displacement u^ the following contributions 
were considered: 

(1) Bending deflection of the yaw leaves as "guided" (i.e. 


having zero slope at the top of the yaw leaves) cantilever beams, 

3 

<S. = (A-7a) 

24ei YWX 


where Iy^ x represents the moment of one yaw leaf about a X axis 
through its centroid. 
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(2) Deflection at orbiter centerline due to axial , 

deformation of the pitch springs. The axial deformations of 
the pitch springs permits a rotation of the orbiter bracket 
about a local X axis. The deflection at the orbiter centerline 


due to this rotation is ' (a+— ) 

6 , = -2 

ApEh 2 


(A-7b) 


where A p is the cross sectional area of one pitch leaf. 

(3) Bending deflection of booster bracket. Each arm of the 
bracket was assumed to act as a cantilever beam with an 


effective length L 2 for bending in the YZ plane* 

3 

■L>, 

6. = 


V + 


6EIBX 


4EI B X 


(A-7c) 


where I nv represents the moment of inertia of one bracket 
arm about an X axis through its centroid. 

(4) The deflection at the orbiter centerline due to end 


slope of booster bracket. 


L 2 2 d dL 2 2 

69 " 4EI^T + 2ET T~ 


(A-7d) 


L BX ‘“ x BX 
The total deflection due to a unit force is 

agg=5 8 +6 7 +6 8 +6 9 . (A-7e) 

The rotation considered due to the unit force is assumed 
to have two components: 


(1) The rotation introduced by the axial deformation of 

(A-8a) 


the pitch leaves, P 

= - 


ApEh 2 

(2) The end rotation of the booster bracket arms bending 
as cantilever beams. 


0 . 


4EI 


BX 


dL 2 

2EI 


(A-8b) 
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The total rotation is 

a 79 = e <* + 0 5 


(A- 8 c) 


The flexibility coefficients a g 7 and 377 are found by 
applying a unit moment abou£ an X axis through the orbiter 
attachment point and' calculating the displacement ug and 
rotation u 7 * By symmetry. 


a 97 = a 79 


(A-9) 


The rotation a-j-j due to the unit moment has two components 
similar to the rotations considered in calculating 379 * Thus 


a 7? = j£_ + L * 


,2 2EI BX * 


(A-10) 


ApEh 1 

Torsion about Z axis: 

The flexibility coefficient a T is the rotation about the Z 
axis due to a unit moment about this axis. Four rotations were 
assumed to contribute to a T : 

(1) The rotation of the pitch leaves about a Z axis through the 
mounting brackets. The pitchleaves are assumed to be clamped 
at their outer ends, 0 


v> 


(A-ll) 


16EI pz 

(2) The rotation introduced at the top of the yaw springs due to 
the yaw springs bending in the XZ plane as cantilever beams. 


*2 = 


2i y 3 


3C EIywy 


(A-llb) 


where C is the distance between the yaw leaves as shown in 
Figure A-l. 



(3) The rotation introduced at the top of the yaw springs due 
to the slope of the arms of the booster bracket bending as 
cantilever beams. 



Using the equations derived above, all of the elements 
in the flexibility matrix for the elastic interface were 
evaluated. The values of the various dimensions used in the 
calculations are tabulated in Table A-I. For the pitch and 
yaw leavesa modulus of elasticity, E = 3 x 10 7 lb/in 2 was used, 
and for the booster bracket, E = 10 7 lb/in 2 . The results of these 
calculations are tabulated in Table A-II. For each nominal 
spring constant the results of two calculations are shown. 

The first calculation in each case was made for effective 
lengths L = 5.5 in. and L 2 = 2.0 in. which were found by 
approximately matching the bracket arm bending deflections to 
the experimentally measured deflection of the arms. The 
second calculation neglects the bending of the bracket arms, 
i.e., the effective lengths are taken as zero. 
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Stiffness Matrix 

The stiffness matrix for the elastic interface in the 
coordinates shown in Figure A- 2 may be written as 


(A-12) 



where all elements outside of the dashed lines are zero. The 
submatrix K z represents the axial stiffness matrix for Z 
displacements; K xz represents the stiffness due to bending in the 

XZ plane; Ky Z represents the bending in the stiffness YZ plane; 
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and ^represents the rotational stiffness about the Z axis. 

Not all of the elements of this matrix are independent. By 
symmetry, of course, Kij = Kji. In addition, equilibrium require- 
ments may also be used to establish other relations between the 
elements in each submatrix. 

Axial and Torsional Stiffness Submatrices: 

The axial an d to rsional problems are similar for analytical 


purposes. For simplicity, only the axial problem will be 
considered; the conclusions reached can be applied to both problems 
The axial stiffness matrix K z has the form 


[K 7 1 = 


K 11 k 12 

_ k 21 k 22 

The four elements are related as will now be shown. Applying the 
definition of K^j as the force at i due to a unit displacement at 
j with all other displacements zero, the elements of Kz are as 
shown in Figure A-3. Then equilibrium of forces, requires 

K 21 = ~ K 11 
k 12 = -k 22 . 

Since, K 21 = -K^, these equations show there is actually only one 

independent element of the axial stiffness matrix. Let this be K z , 

K Z -K z 


then in general. 


[K z 3 - 


-Kr 


Kr 


Similarly, the torsional stiffness matrix is 


[Kip J - 


Kr 


— Kn 


-K»p 

Kip 


(A-13) 


(A-14) 


Bending Stiffness Matrices: 

The bending stiffness submatrices [K XZ J and [Ky Z ] are 
similar and only one matrix will be considered. The matrix 
IK XZ J is defined, from equation (A-12) , as 
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t R xzl 


k 33 

k 34 

*35 

r 36 

k 43 

k 44 

*45 

r 46 

k 53 

k 54 

k 55 

k 56 

k 63 

*64 

*65 

k 66 


(A- 15) 


Because of symmetry there are 10 unknown values in the array. 
When equilibrium requirements are considered it turns out that 
only three of the values are independent. This follows from 
applying the definition of K^j as shown in Figure A-4. 
Equilibrium of forces and moments gives the following eight 
relations , 

\ Forces =0 
\ Moments =0 

0 

(Figure A4-a) 


k 53 

+ k 63 = 0 


K 33 + 

k 43 + lk 53 = 

= 0 

K 54 

■f K 64 - 0 


K 34 + 

K 44 + lk 54 

= 0 

K 55 

+ K 65 - 0 


K 35 + 

k 45 + lk 55 

= 0 


(Figure A4-b) 


(A- 16) 


(Figure A4-c) 


k 56 + K 66 = 0 
k 36 + k 46 +.1*56 = 0 


(Figure A4-d) 


With these eight relations and taking symmetry into account, 
the matrix above in equation (A-14) can be written as 

V 40- . /VJfc':- 



(A-17) 


K 33 

K 34 

_ K 33 + K 34 

k 33 + k 34 



L 

L 

K 34 

K 44 

K 34 + K 44 

k 34 + k 44 



L 

L 

k 33 + k 34 

k 34 + K 44 

K 33 + 2k 34 A K 44 - 

. k 33 + 2k 34 + k 44 

L 

L 

L 5 

Li 

k 33 + k 34 

K 34 + k 44 

k 33 + 2k 34 + k 44 

k 33 + 2k 34 + k 44 


L L L 5 L 5 

The matrix [Kygl may be found in the same form by changing 
subscripts appropriately. 

As a result of symmetry and equilibrium requirements it 
has been shown that there are only 8 independent stiffness 
coefficients: K z and K»p and 3 elements each from the submatrices 

Kxfc and Ky Z . 

Relations between the Stiffness and Flexibility Coefficients 
To determine the stiffness coefficients in terms of the 
flexibility coefficients the definition of a^j is used. For the 
interface in extension and torsion this means simply that 

K Z = (A- 18) 

and 

1 

k T = a T . (A-19) 

For the interface in bending this means that one end is rigidly 
fixed and a unit force and a unit moment are applied at the free 
end. The flexibility coefficients can then be used to write, 
assuming symmetry. 
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a 33 a 35 


a 35 a 55 


where P's denote the applied forces. In terms of the stiffness 
matrix of the interface in bending 

fPoi .k, c i r 


K 33 K 35 

K 35 K 55 


Thus if the flexibility coefficients are known, 
F K 33 K 35*1 f a 33 a 35] 1 


K 35 K 55 a 35 a 55 


a 33 a 55 " a 35 2 


55 

” a 35 

35 

a 33 


hence , 


‘33 = 


a 33 a 55 " a 35 2 


_ a 35 

K 35 - a 33 a 55 “ a 3 5 2 


(A- 20) 


55 a 33 a 55 “ a 35 2 

Using these relations and the equations of equilibrium, (A-16) 
found above, the stiffness matrix [K xZ ] can be found in terms 


of the flexibilities by: 
f a cc 


(A-21) 


La 35” a 55 


La 35“ a 55 a 55“ 2La 35 +L2a 33 a 35” La 33 “’ a 35 +La 33 


a 33 

" a 33 


” a 33 

a 33 


a 33 a 55 _a 35 2 " a35 a 35” La 33 a 33 ” a 33 

a 35 " a 35 +La 33 " a 33 a 33 

• - -J • 

A similar relation for [ K yz ] can f° un <3 by appropriately changing 
the subscripts. 
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Equations (A-18), (A-19) , and (A-21) were used to calculate 
the stiffness coefficients corresponding to the flexibility 
coefficients given in Table A-II. These stiffness coefficients 
are tabulated in Table A-III. As in Table A-II, for each 
nominal spring constant the stiffness coefficients were 

ca leu la ted “for - two cases of~~ef f ec ti ve~l~eng ths o f the— boos ter — 

bracket arms . 

In this appendix, equations have been derived which may be 
used to calculate the flexibility matrix for the elastic 
interface. In addition explicit relations have been derived 
which may be used to calculate the elements of the corresponding 
interface stiffness matrix directly from the elements of the 
flexibility matrix. Numerical results for the flexibility and 
stiffness matrices for each nominal case of the elastic 
interface were presented. 
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APPENDIX B 

Rigid Body Vibration Analysis 

In this appendix a description of the rigid body vibration 
analysis will be given. Two analytical models were used for 

this— analysis-. — -I-n— the— f-i-rst— analysis-,— the— motion— was— assumed 

to be restrained to the pitch plane with the longitudinal 
displacement' prohibited. This gave a simple four degree of freedom 
model. In the second analysis, general three dimensional motion 
was permitted. The booster and orbiter each were allowed to have 
three translations and three rotations so that the system possessed 
twelve degrees of freedom. 

Plane Motion 

The mathematical model for the case of plane motion is shown 
in Figure B-l. The generalized coordinates are the vertical 
translations Z l and Z 2 of the mass center of each component 
and the rotations 0 X and 0 2 about the mass center. Each elastic 
interface is represented as a scalar spring having a single spring 
constant. 

For free vibrations , the equations of motion for small 
oscillations takes the form 

[M] (q) + [K] {q} = 0 (B-l) 

where [M] denotes the mass matrix and [K] denotes the stiffness 
matrix. For the present problem the generalized coordinates 
(q } are given by ; 

{q} = |_Z, Z, 0, jT (B _ 2 ) 
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The mass and stiffness matrices were obtained by writing 
the kinetic and potential energies of the system. Since the 
generalized displacements are measured from the mass center of 
each body the mass matrix is diagonal and was found to be 



(B-3) 


From the geometry shown in Figure B-l, the relative 
displacement of each spring was written in terms of the 
generalized coordinates. With these displacements, the 
potential energy of the system was derived and written in 
matrix form to yield the stiffness matrix. For the present 
problem these calculations gave 

(B-4) ' 


[K] = 


-(K-j+Kj) 


- (K x +K 2 ) (Kj+I^) 


X K, -XK 0 -X K, +X K 0 
alc2 a 1 c 2 


-X b K i+ X d K 2 x b K l~ X d K 2 


X a K l- X c K 2 


-X a K 1+ X c K 2 


- X b K l +X d K 2 


x b x l-X d K 2 


X a 2K l +X c 2K 2 -< X a X b K l +X c X d K 2 ) 
-(Wl + X c X d K 2> WW 


Using the mass and stiffness matrices derived above the 
eigenvalue problem 


[K] {u} = u) 2 [M] {u} (B-5) 

was then solved to give the system's natural frequencies w and 
mode shapes {u}. Since two "free-free" modes are possible, 
two zero frequencies were found as well as two frequencies 
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corresponding to the system's modes of . vibration . 

General Motion 

Figure B-2 shows the mathematical model for the general 
case of three dimensional motion. Each component was permitted 
to have six degrees of freedom: three translations of each mass 
center and three rotations about each mass center. The positive 
conventions for these coordinates are shown in Figure B-2. 

Each spring interface was represented by a twelve by twelve 
stiffness matrix. 

The generalized coordinates in this case may be partitioned 
into orbiter and booster displacement sets as 
(q) = 


u 0 

U B 


(B-6) 


where the displacements of the orbiter are 


T 

9 


(u 0 > « [u 10 u 20 U 30 e 10 e 20 ®30j • (B-7a) 

and the displacements of the booster are 

T 

{u B } = L UlB U2b U3b ® 1b 028 03B J ’ (B-7b) 

For the fuselages, the coordinate axes through the mass center 
are principal axes because of symmetry. With the wings attached, 
this is also a good approximation since the wing mass is relatively 

small. Thus the displacement coordinates will be assumed to be aligned 
and referenced to the principal axes. Then the mass matrix is 


[M] = 


”0 


1 


.“B 


(B-8a) 


where the diagonal elements of [Mq] are 
^0 ^0 M Z0 I X0 


*Y0 I Z0 ' 


(B-8b) 
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are 


and the diagonal elements of [Mgl 

M XB M ZB r XB X YB I ZB * (B-8c) 

All off-diagonal elements are zero. M x q denotes the mass 
associated with the X displacement of the orbiter; I^o denotes 
the mass moment of .inertia of the orbiter about an X axis through 
the center of mass. The other terms associated with the Y and Z 
subscripts may be similarly interpreted. The analogous quantities 
for the booster are denoted with the B subscript. 

The stiffness matrix for the system is the matrix sum of 
the stiffness matrices for each interface where each matrix 
is a twelve by twelve. Thus 

[K] = [Kj ] + [K 2 ] (B-9) 

where [K x ] represents the interface connecting points a and 
b, and [K 2 ] represents the interface connecting points c and d 
in Figure B -2 . For simplicity, the derivation of only the [K ] 
matrix will be given. It will be indicated later that the 
elements of the [K 2 ] matrix may be obtained from [Kj] by making 
appropriate substitutions. 

To derive the stiffness matrix [Kj] it is desirable to begin 
with the set of coordinates used in Appendix A to represent 
the interface. These coordinates, and the associated stiffness 
matrix for the interface in these coordinates are shown in 
Figure B-3. The sub-matrices in Figure B-3 may be identified 
as follows: 
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[K z J represents the axial stiffness of the interface to 
relative displacements in the z direction, two by two; 

[K xz ] represents the stiffness of the interface to relative 
displacements and rotations in the xz plane, four by four; 

[K ] represents the stiffness of the interface to relative 
displacements and rotations in the yz plane, four by four; and, 
[K T ] represents the torsional stiffness of the interface to 
relative rotations about the z axis, two by two. 

It is desired to represent the stiffness of the interface in 

terms of the displacements at the center of mass. This will be 
done by writing the potential energy of the interface first in 
terms of the coordinates of points a and b, and then making a 
coordinate transformation. The potential energy may be written 


as 


where 


1 

u a 

T 

K aa 

K ab 


u a 

2 

u b 


K ba 

K bb 


u b 


( u a ) ~ [ u la u 2a u 3a 0 la 0 2a 0 3a J 7 


(B-10) 


(B-lla) 


and 

(u b > = |u lb u 2b u 3b 0 lb 0 2b © 3b J f (B-llb) 

The submatrices in equation (B-10) are found by identifying the 
(ua) and {%) displacement components with the displacement 

components shown in Figure B-3, and then appropriately reordering 
rows and columns in the stiffness matrix. The elements of the 
submatrices are defined in Figure B-4 where K^j denotes an element 
in the matrix [K X z3 and K'^j denotes an element in [Kyz3 • 
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Refering to Figure B-2 the transformation equations 
between the coordinates at points a and b and the generalized 
coordinates may be written. Assuming small rotations 

{u a } = [.8 a ] {uq} (B-12a) 

where 


[B a ] = 


Similarly, 


where 


[BvJ = 


0 

1 


0- 

0 


0 0 1 
0 0 0 


- 0 - 

Zc 

0 

1 


— 2 . 

c 

0 
X = 


0 

_X c 

0 

0 


0 

0 

0 

1 0 

1 

0 

0 

0 

0 

1 

1 0 

0 

1 

{u b } 

= C6 b ] 

{u B } 




1 

0 

0 

1 

1 .0 

Zb 

0 


0 

1 

0 

J -Zb 

0 

-Xb 


0 

0 

1 

! 0 

Xb 

0 


0 

0 

0 

! 1 

0 

0 


0 

0 

0 

! 0 

1 

0 


_ 0 

0 

0 

l 0 

0 

1 


(B-13) 


(B-14) 


(B-15) 


Substituting equations (B-12) and (B-13) into equation 

it 

(B— 10 ) gives 


Vl = 1 


_ u o 

T 

£■• id 
cn 

i 

! o 


K aa 

! K ab " 


"B a 

0 


u 0 

U B 


0 

i . 


" K ba 

T K bb". 


0 

’®’b’ 


. U B . 
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From this, by inspection'. 


[KJ 




T 

*b 


K 


aa 

K ba 


K 


ab 


K 


bb 


0 

e b 


Finally, after performing the multiplications. 


[K x ] = 


"g T 

Pa 

K aa 

e a 

1 ft T 
1 “a 
-1 

K ab 

^b 

T 

_*b 

Kba 

S a 

! *1 

Kbb 

&b 


(B-16) 


Each of the submatrices is of order six by six. When these 
operations were carried out the final matrix was obtained. The 
non-zero elements of the upper triangle of this matrix are 
listed in Figure B-5. The elements of the lower triangle may 
be found by symmetry. 

The elements of the [K 2 ] stiffness matrix may be found 
from the [K x ] matrix by replacing: 

X a by -X c 

z a hy z c 
X a by -X d 

z b ky z d • 

These replacements were found by writing coordinate transformations 
similar to equations B-12 and B-14 for points c and d. 

With the mass matrix and the stiffness matrix for the 
system as given in equations B-8 and Figure B-5, the free 
vibration modes were found by solving an eigenvalue problem 
of the type given in equation (B-5). In this case there are 
six "free-free" modes and six corresponding zero frequencies 
in addition to six natural modes of vibration. 
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The mass and inertia properties for the plane motion 
analysis and the general motion analysis are tabulated in 
Table BI. These values were found using the corresponding 
NASTRAN finite element model in each case. The rigid body 
mass and inertia properties are found in NASTRAN by requesting 
execution of the Grid Point Weight Generator. The Grid Point 
Weight Generator is requested by use of the bulk data PARAM 
card. 
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Station 15.7 36.0 73.2 96.7 117.4 150.0 

(inches) 


30.7 54.4 81.7 102.7 139.9 


(b) Elevation view. 


Figure 2.- Schematic of 1/15 scale space shuttle model 







Fuselage 



Figure 4 . - Schematic of fuselage-wing interface 





Figure 5.- Schematic of booster delta wing. 








Figure 8.- Schematic of booster fuselage-wing interface geometry. 
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^^■^ure 10- — Schematic of booster airplane idealization. 




Frequency, HZ 


Mode 

Liftoff 

Max Q 

Burnout 

~ 1" 

2 

3 

' 36.8 
117.2 
198.7 

53.7 

137.9 

252.3 

83.9 

243.7 

464.1 





Figure 12.- Bending inodes of booster fuselage. 
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(a) Plan view of delta wing. 
Z 




- X 

• (c) Undeformed wing as seen by plotter. 

Figure 14.- Plotter orientation for delta wing mode shape plots 





(b) Mode 2, 26.4 HZ. 



Figure 15.- Mode shapes of booster delta wing. 
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Figure 15 (concluded).- Mode shapes of booster delta wing. 
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I 


(b) Mode 2, 29.3 HZ. 



(c) Mode 3, 31.0 HZ. 


Figure 16.- Mode shapes of orbiter delta wing. : 
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(a) Analytical model 







Figure 17.- Mode shapes of elastically connected fuselages. 

Liftoff weight distribution. 

10 2 lb/in nominal springs. 
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;p] 



r 

Rigid body 
NASTRAN 


j J 


Mode 2^/ 

s ' flode 1 



Nominal spring constant, lb/m. 


Figure 19.- Variation of NASTRAN and rigid body analyses 
frequencies with nominal spring constant. 



■<v 




Frequencies 

Rigid body - 3.44 HZ 
NASTRAN - 3.53 HZ 


(a) First mode. 



Rigid Body 
NASTRAN 


Figure 20.- Comparison of NASTRAN and rigid body analysis modes. 

10 2 lb/in nominal spring . 
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(b) Displacement components for mode shapes. 


Figure 21.- Plotter view of elastically connected fuselages 


(a) Mode 1, 11.2 HZ 


T3=- . 4 


T2=. 8 


Tl=-. 1 


T2=-l. 


(b) Mode 2, 17.3 HZ. 


Tl=- . 1 
T3=l. 


Cc) Mode 3 , 24.9 HZ 


Piqure 22.- Mode shapes of elastically connected fuselages 
Theoretical flexibility matrix. Nominal 10 lb/xn. 















(b) Plotter view of undeformed model. 

Figure 25.- Plotter orientation for complete model mode shape plots. 
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Figure 26 (concluded).- Mode shapes of booster airplane. 
Liftoff weight distribution. 
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rH (N H 












(f) Mode 6, 30.9 HZ. 

Figure 28 (continued) .-Mode shapes of complete model. 
10 3 lb/in theoretical flexibility matrix. 
Liftoff weight distribution. : 
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T3=-l. 

(h) Mode 13, 65.6 HZ. 


Figure 28 (continued).- Mode shapes of complete model. 
10 3 lb/in theoretical flexibility matrix. 
Liftoff weight distribution. 
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(j) Mode 18, 


78.7 HZ. 



Figure 28 (concluded) Mode shapes of complete model. 
10 3 lb/in theoretical flexibility matrix. 
Liftoff weight distribution. 
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Figure B-l.- Plane motion rigid body model. 



(a) Orbiter-booster system 


u 30 



U 3B 



IB 


(b) Orbiter displacement (c) Booster displacements. 


Figure B-2.- General motion rigid body model. 
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*Kjl j denotes an element of [K XZ J • 
** K 'ij denotes an element of [Ky Z ], 


Figure B-4.- Stiffness matrix [K^] in attachment coordinates. 
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*, 1 - *33 

33 

* 1 ' 1 =K 23 + ^ K 3 ‘ 

* 1,11 2 

*2 2 =K ’ 33 

; K ', , + Z a K 33 
*2,4 13 

- -X * ' 33 
* 2 . 6 " a 

* 2,8 =* 3 * k . . 
= k< 23 * ZfaK 34 

* 2,10 

« -X b *' 34 
*2.12 * 

*3,3 = * A 


*3,5 


« X a *a 


*3,9 = 

= -X b *A 

* 3 ' , » K -33 + i *» K ‘ t3 +K 

^■;i a ^ 33-^ 33 
K 4 ’,B= K '3* + ' aK ' 13 


l B a ***' 34 +K ' 33 

**,«“ *• 23 


- & b *' 14 


* *’ 


figure 


B 


^ 5> - -Elements 


of l*ll in 


genera 


lized 


coordinates 
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K 4 , 12 = -Z a X b K *3 4 " X b K 'l4 
k 5 , 5 = Z a 2 K 3 3 + x a 2 K Z ’ 2Z a K 13 + K 11 
K 5,7 K 1 4 " Z a K 34 
K 5 , 9 - -X a K 2 

K 5 , 11 = -Z a K 2 3 - Z a Z b K 34 ~ X a X b K Z + K 12 + Z b K 14 

k 6 , 6 = x a 2 K 33 + k T 

K 6 , 8 = -X a K *34 

k 6 , 10 = -x a K *2 3 + x a z b K< 34 

k 6 , 12 = x a x b K ' 3 4 “ K T 
k 7 , 7 = k 44 

k 7 , 11 = k 24 + z b k 44 
k 8,8 = K *44 

k 8 , 10 =k '24 “ z b k *44 
k 8 , 12 = ~ x b K *44 
k 9 , 9 = K z 

K 10,10 = z b 2K '44 " zz b K 24 +K '22 
k 10,12 = x b z b K 44 ~ x b K '24 
X 11 , 11 ~ ^ z b x 24 + x b* X Z + x 22 + ^b x 44 
k 12,12 = X b 2 K 44 + k T 


Figurs B - 5 (continued) • ** Elements of [K^] in generalized coordinates 
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TABLE I 

Booster and orbiter simulated propellant weights 


(a) Booster add-on weights 
(Pounds) 


Weight No. 
(See Figure 2) 

Liftoff 

Max Q 

Bif?nout 

1 

138.7 

138.7 

0 

2 

22 . 4 

22.4 

0 

3 

80.6 

0 

0 

4 

13.9 

0 

0 


(b) Orbiter add-on weights 
( Pounds ) 


Position 

Liftoff 

Empty 

Forward 

60.4 

0 

Aft 

9.9 

0 
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Table I I 

Booster fuselage 
Element properties 


Grid point (typical) 
12 3 


TO" 




Cross Section Properties 

NSM ( lb-sec 2 /in/in) 

Element 

A ( in 2 ) 

I(in 4 ) 

J(in 4 ) 

Liftoff 

Max Q 

1 

1. 305 

0.7882 

1.576 

0 

0 

2 

1.261 

3.964 

7.928 

2.243-2* 

0 

3 

1.707 

9.843 

19.69 

0.965-2 

0.965-2 

4 

2.072 

17.61 

35.22 

0 

0 

5 

2.092 

18.21 

36.42 

0 

0 

6 

2.113 

18.77 

37.54 

0 

0 

7 

2.203 

21.27 

42.54 

0 

0 

8 

2.416 

28.04 

56.08 

0 

0 

9 

2.602 

35.04 

70.08 

0 

0 


*The notation 2.243-2 represents 2.243 x 10 
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Table III 
Booster fuselage 
Concentrated mass properties 


Liftoff 



Nomina' 

7 

. Values 

Test Values 

Grid 

Station 

M 

1 V V 

- 

M 

-*-xx 

Point 

x ( in) 

lb-sec 2 /in 

in-lb-sec 2 

lb-sec 2 /in 

in-lb-sec 2 

1 

■ ■■■bh 

0.233-3 

0.0023 

0.233-3 

0.0023 

2 


0.207-2 

1.366 

0.207-2 

1.47 

3 

• 

0.184-2 

4.937 

0.174-1 

4.28 

4 

| | | | 

0.298-2 

3.636 

0.139-1 

2.86 

5 

66.0 

0.360-1 

0.906 

0.360-1 

0.83 

6 

80.97 

0.362-2 

0.135 

0.362-2 

0.135 


87.0 

0.373-2 

0.148 

0.373-2 

0.148 


101.7 

0.580-1 

1.734 

0.580-1 

1.69 


124.2 

0.455-2 

0.325 

0.455-2 

0.325 

■■■Hi 

134.8 

0.142-1 

0.345 

0.246-1 

0.345 




Nominal 



Max Q 

Burnout 

Grid 

Station 

M 

X xx 

M 

J xx 

Point 

x (in) 

lb-sec 2 /in 

rn-lb-sec 

lb-sec 2 /in 

in-lb-sec 2 

1 

0 

0.233-3 

0.0023 

0.233-3 

0.0023 

2 

11.0 

0.207-2 

0.0172 

0.207-2 

0.0172 

3 

20.3 

0.184-2 

3.588 

0.184-2 

0.110 

4 

57.5 

0.298-2 

3.636 

0.298-2 

0.158 

5 

66.0 

0 

0.107 

0 

0.107 

6 

80.97 

0.362-2 

0.135 

0.362-2 

0.135 


87.0 

0.373-2 

0.148 

0.373-2 

0.148 


101.7 

0.580-1 

1.734 

0 

0.240 


124.2 

0.455-2 

0.325 

0.455-2 

0.325 


134.3 

0.142-1 

_ 

0.345 

0 .142-; 

0.345 
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Table IV 
Orbiter fuselage 
Element properties 


Grid point 
20 21 


(typical) 

22 



24 


25 


26 27 


a 

© * 

~© 

— 

'(0>t ® ■ 


TW> 



V X 



Elemen 

t numbe: 

c (typical) 


Element 

Cross section properties 

NSM ( lb-sec 2 /in/in) 

-/V'V' . jfl \ 

I(in 4 ) 

J (in 4 ) 

Liftoff 

Empty 

hh 

1.169 

3.160 

6.32 

0 

0 


1.458 

6.138 

12.28 

1.042-2 

0 


1.272 

4.072 

8.14 

0 

0 

1 

1.272 

4.072 

8.14 

0 

0 


1.272 

4.072 

8.14 

0.214-2 

0 

mi 

1.272 

4.072 

8.14 

0 

0 

mam 

1.272 

4.072 

8.14 

0 

0 
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Table VI 

Booster wing grid point coordinates 


Coordinate System 1 


Grid Point 

X i ( inches) 

Y (inches) 

Z j ( inches) 

1 

41.33 

0 

o 

2 

20.66 

0 

0 

3 

0 

0 

0 

4 

41.33 

-9.42 

1.62 

5 • 

20.66 

-9 . 42 

1.62 

6 

0 

-9.42 

1.62 

7 

41.33 

9.42 

1.62 

8 

20.66 

9.42 

1.62 

9 

0 

9.42 

1.62 

34* 

41.33 

-5.82 

1.00 

35* 

20.66 

-5.82 

1.00 

36* 

0 

-5.82 

1.00 

37* 

41.33 

5.82 

1.00 

38* 

20.66 

5.82 

1.00 

39* 

0 

5.82 

1.00 


♦Grid points 34-39 were used only for attaching the 
wing to the fuselage. 


Coordinate System 2 


Grid Point 

X 2 { inches) 

Y 2 (inches) 

Z 2 (inches) 

10 

10.0 

41.33 

0 

11 

10.0 

25.82 

0 

12 

10.0 

10.31 

0 

16 

20.0 

41.33 

0 

17 

20.0 

30.97 

0 

18 

20.0 

20.62 

0, 

22 

30.0 

41.33 

O' 

23 

30.0 

36.13 

0 

24 

30.0 

30.94 

0 

28 

35.0 

41.33 

0 

29 

35.0 

38.71 

0 

30 

35.0 

36.10 

0 
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Table VI continued 


Coordinate System 3 


Grid Point 

X 3 ( inches)' 

Y' s (inches) 

Z (inches) 

13 

41.33 

10.0 

0 

14 

25.82 

10.0 

0 

15 

10.31 

10.0 

0 

19 

41.33 

20.0 

0 

20 

30.97 

20.0 

0 

21 • 

20.62 

20.0 

0 

25 

41.33 

30.0 

0 

26 

36.73 

30.0 

0 

27 

30.94 

30.0 

0 

31 

41.33 

35.0 

0 

32 

38.71 

35.0 

0 

33 

36.10 

35.0 

0 
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Table VII 

Booster wing element properties 



A 

11 

12 

J 

112 

Element 

in 2 

in 4 

in 4 

in 4 

in 4 

1,11,24,10,20,33 

0.307 

0.596-2 

0.129-1 

. 

0.106-1 

0 

2,12,25,9,19,32 

0.419 

0.299-1 

0.169-1 

0.119-1 

0 

3,13,26,8,18,31 

0.559 

0.118 

0.200-1 

0.136-1 

0 

4,14,27,7,17,30 

0.695 

0.299 

0.219-1 

0.152-1 

0 

40,41 

0.714-1 

0.281-3 

0.720-2 

0.900-4 

-0.895-3 

50,51 

0- 714-1 

0.396-2 

0.720-2 

0.900-4 

+0.895-3 

42,43 

0.936-1 

0.296-2 

0.980-2 

0.111-3 

-0.308-2 

52,53 

0.936-1 

0.296-2 

0.980-2 

0.111-3 

+0.308-2 

44,45 

0.116 

0.114-1 

0.114-1 

0.139-3 

-0.622-2 

54,55 

0.116 

0.114-1 

0.114-1 

0.139-3 

+0.622-2 

5,15,28,6,16,29 

1. 33 

0.857 

0.143 

0.320-1 

0 

35,36,37,38 

0.25 

0.130-2 

0.208-1 

0.833-1 

0 
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Table VIII 

Orbiter wing grid point coordinates 


Grid point 

X (inches) 

Y (inches) 

Z (inches) 

1 

2 8.36 

0 

0 

2 

14.36 

0 

0 

3 

0 

0 

0 

4 

28.36 

4.928 

0.8468 

5 

14.36 

4.928 

0.8468 

6 

0 

4.928 

0.8468 

7 

28.36 

-4.928 

0.8468 

8 

14.36 

-4.928 

O'. 8468 

9 

0 

-4.928 

0.8468 

10 ' 

28.36 

14.78 

2.540 

11 

19.34 

14.78 

2.540 

12 

10.31 

14.78 

2.540 

13 

28.36 

-14.78 

2.540 

14 

19.34 

-14.78 

2.540 

15 

10.31 

-14.78 

2.540 

16 

28.36 

24.64 

4.234 

17 

24.49 

24.64 

4.234 

18 

20.63 

24.64 

4.234 

19 

28.36 

-24.64 

4.234 

20 

24.49 

-24.64 

4.234 

21 

20.63 

-24.64 

4.234 

22 

28.36 

29.57 

5.081 

23 

27.08 

29.57 

5.081 

24 

25.78 

29.57 

5.081 

25 

28.36 

-29.57 

5.081 

26 

27.08 

-29.57 

5.081 

27 

25.78 

-29.57 

5.081 

44* 

28.36 

-3.94 

0.676 

45* 

14.18 

-3.94 

0.676 

46* 

0 

-3.94 

0.676 

47* 

28.36 

3.94 

0.676 

48* 

14.18 

3.94 

0.676 

49* 

0 

3.94 

0.676 


♦Grid points 44-49 were used only for attaching 
wing to the fuselage. 














ies 



ment properti 


12 ^ 

J 

112 

in 14 

in 4 

in 4 

0.160 

0.293-2 

0 

0.998-2 

0.352-2 

0 

— 0 7751-2 

0.302-2 

o 

0.509-2 

0.326-2 

0 

0.260-2 

0.260-2 

0 

0.848-2 

0.109-3 

• 

-0.618-3 

0.848-2 

0.109-3 

-0.618-3 


Table X 


Theoretical and Experimental Flexibility 
Matrices for elastic interface. 


Nominal 10 3 Case 


Eq. 

IJASTRAIJ 

7 . . * 

_ z ij 

(AI-l ) 

Theoretical 

Experimental 

a 55 

Z 11 

0.1041-2 

0.1035-2 

a 53 

z 15 

0.1344-3 

0.1324-3 

a Q0 

z 22 

0.1055-2 

0.1215-2 

a 97 

Z 24 

-0.0775-4 

-0.216-4 

a ll ~ a z 

Z 33 

0.9962-3 

0.976-3 

a 77 

CSJ 

0.1112-5 

0.233-5 

a 33 

Z 55 

0.1814-4 

0.1807-4 

a T 

Z 66 

0.5263-5 

0.721-5 


all other elements are zero 












Table XI Resequencing of grid points for complete model. 
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Table XII Natural frequencies of booster and orbiter fuselages, 
(a) Booster fuselage frequencies, HZ. 




Weight distribution 



Liftoff 

Max Q 

Burnout 


Mode 

Theory 

Exper. 

Theory 

Theory 

Exper . 


1 

36.8 

36. 3 a , 36. 8^ 

53.6 

83.9 

73.7 

Bending 

2 

117.2 

95. 2 a , 98. 8 b 

137.9 

242.7 



3 

198.7 

181.9, 177.3 

252.3 



Torsion 

1 

114.3 

88.9 

116.8 


mi 


2 

243.4 

161.8 

247.5 



Axial 

1 

240.2 


298.9 



a Pitch excitation 
b Yaw excitation 


(b) Orbiter fuselage frequencies, HZ. 



Weight distribution 

Liftoff 

Empty 


Mode 

Theory 

Exper. 

Theory 


1 

103.5 

103. 4 a , 102. 7 b 

■B— 

Bending 

2 

229.7 

205°, 201 c 



3 

460.3 


mi 


a Pitch excitation 
k Yaw excitation 
Approximate 
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Table XIII Natural frequencies of free-free booster and orbiter v/ings 



Booster wing 

Orbiter wing 

Mode 

Theory 

Exper . 

Theory 

Exper. 

1 

24. 

6 

'26.4 

26.3 


2 

26. 

4 

24.8 

29.3 


3 

27. 

3 

26.4 

31.0 

28.5,32.2 

4 

28. 

0 

26.0 

61.3 


5 

49 . 

4 


63.2 

66.9 

6 

54. 

4 

60.4 

67.5 

74.5 

7 

65. 

0 


69.1 


8 

69. 

4 


73.4 


9 

69. 

6 

73.8 


100.7 

10 

96. 

1 



112.5 

11 

101. 

0 


118.4 


12 

107. 

3 



120.5 

13 

113. 

2 


12 5.7 


14 

125. 

2 


128.9 


15 

142. 

5 


133.4 


16 

142. 

8 


139.3 


17 

152. 

2 


177.1 


18 

155. 

2 


183.8 


19 

182. 

9 


190.4 


20 

194. 

2 


223.2 
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Table XVI Comparison of frequencies for elastically connected fuselages. 

Nominal 10 3 lb/in interface. 

Liftoff weight distribution. 


Nominal Interface 
Scalar springs 3 

Interface Flexibility Matrix 

Theoretical 'Matrix* 3 

Experimental Matrix* 3 

HZ 

Motion 

HZ 

Motion 

HZ 

Motion 

10.7 

Pitch 


Pitch 


Pitch 

10.7 

Yaw 


Yaw 


Yaw 

25.9 

Pitch 


Pitch 


Pitch 

25.9 

Yaw 


Yaw 


Yaw 

36.3 

Pitch 

37.4 

Pitch 

37.3 

Pitch 

36.3 

Yaw 

40.1 

Yaw 

38.5 

Yaw 

93.8 

Pitch 

75.2 

Pitch-axial 

64.3 

Pitch-axial 

93.8 

Yaw 

96.2 

Yaw 

88.4 

Roll 

110.5 

Pitch 

104.9 

Pitch 

95.5 

Yaw 

110.5 

Yaw 

110.4 

Roll-yaw 

102.6 

Pitch 

191.0 

Pitch 

113.1 

Yaw 

112.9 

Pitch 

191.0 

Yaw 

118.7 

Pitch 

116.7 

Pitch 



131.9 

Roll 

i23 . 9 

Roll 


^Test weights. Inverse power method with consistent mass matrix. 
to Test weights. Givens method with lumped mass matrix. 















Table XVII Frequencies of booster and orbiter airplanes. 

(a) Booster airplane. 


Frequencies, HZ 



Predominate 

motion 


Coupled fuselag 

e and wing modes 

Wine modes 

Pitch 9 

Yaw b 

Symmetric 0 

Antisymmetric 0 

30. 0 d 
3 3.9' 
122.3 
128.1 

30. 7 d 

TIM 

42.53 
63. 8 d 

69.5 

67.5 
119.8 


5* Pitch refers to displacements in the XZ plane. 
b Yaw refers to displacements in the YZ plane. 

° Symmetry and antisymmetry refer to the XZ plane. 
d The mode shapes for these frequencies are plotted in Figure 28. 


(b) Orbiter airplane. 


Frequencies, HZ 



Predominate 

motion 


Coupled fuselage and wing modes 

Wing modes 

Pitch 9 

Yaw b 

Symmetric 0 

Antisymmetric 0 

79. 6 d 

76. l d 

3 3. l d 

41. 6 d 



70. 6 d 

71. 7 d 



132.7 

128.1 



138.8 

141.7 



152.2 

149.7 




198.9 


a Pitch refers to displacement in the XZ plane. 
b Yaw refers to displacements in the YZ plane. 
c Symmetry and antisymmetry refer to the XZ plane. 
d The mode shapes for these frequencies are plotted in Figure 28. 
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Table XVIII Frequencies of complete model 
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Table A-I . 

Spring Interface Dimensions 


Nominal Spring Constant 
(lb/in ) 

t inches 

ty inches 

10 2 
10 3 
10 14 
10 s 

0.0188 
0.0406 
0.0874 
0 ..1.8 8.2 

.0237 
.0511 
. 1101 
...2 37-1- 

' 




Dimension 

Value 

inches 

a 

6.5 

h 

1.6 

ip 

2 . 0 a 

ly 

2.0 b 

Lx 

5.5 

L 2 

2.0 

d 

9.0 

C 

8.0 


a In equations A-3c, 4b, 6, 7b, 8a, 11 an effective 
length £ p = 3.0 inches was used. 

b In equation A-llb and 11c an effective length Zy = 2.75 
inches was used. 
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SPRING INTERFACE FLEXIBILITY MATRIX 
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